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Abstract 

A two-loop calculation of the renormalization group /^-function in a momen- 
tum subtraction scheme with massive quarks is presented using the background 
field formalism. The results have been obtained by using a set of new generalized 
recurrence relations proposed recently by one of the authors (O.V.T.). The be- 
havior of the mass dependent effective coupling constant is investigated in detail. 
Compact analytic results are presented. 
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1 Introduction 



Particle masses are amongst the most important physical parameters and in many cases 
their meaning and definition by thresholds (e.g. lepton masses), symmetry breaking 
parameters (current quark masses, neutrino masses) or scale parameters is quite clear. 
For particles which exist as free or quasi-free states a definition by the pole mass is 
most natural and has an unambiguous meaning. The definition of quark masses, in 
particular for the light quarks, allows for a lot of freedom, mainly because the pole mass 
is not directly observable due to the confinement property of QCD. Nevertheless, quark 
masses play a crucial role for the effective behavior of strong interactions at a given 
scale. The purpose of the present calculation is a precise understanding of the quark 
mass dependence of QCD, more specifically, of the effective coupling constant a s (Q 2 ) = 
<7g/(47r), the most important quantity in the description of strong interactions. These 
considerations are important for a better understanding of the decoupling of heavy 
particles and of the relationship between QCD with massive quarks and QCD in the 
MS scheme where effective theories with different number of (light) flavors must 
be matched at the different quark thresholds. 

When the on-shell renormalization scheme is not adequate we either may use a minimal 
subtraction scheme (MS or MS ) or some version of a momentum subtraction scheme 
(MOM ), defined by the condition that the radiative corrections of an appropriate set 
of quantities vanish at a certain (off-shell) momentum configuration. While the MS 
[|J scheme is technically simple and respects the Slavnov- Taylor identities the MOM 
scheme is more physical since it respects the decoupling theorem A serious short 
coming of the standard MOM schemes 0, however, is the fact that they spoil the 
validity of the canonical form of the Slavnov- Taylor identities. An elegant way out 
of this difficulty is the use of the so called background field method (BFM) 0. The 
latter takes advantage of the freedom to chose a gauge fixing function in a particular 
way, namely, such that the canonical Slavnov- Taylor identities remain valid also after 
momentum subtractions. The gauge invariant physical quantities are not affected by 
the gauge fixing, however, the "background field gauge" selects a particular representa- 
tive of the gauge variant off-shell amplitudes. The restauration of the Slavnov- Taylor 
identities in the BFM is achieved solely by changing the vertices with external gluons 
appropriately. For further details and for the Feynman rules of QCD in the background 
field (BF) formalism we refer to 

In Ref. the renormalization group (RG) /3-function of QCD was evaluated at one- 
loop order in the background field approach using the MOM scheme. In the present 
article we extend this analysis to a complete mass dependent two-loop calculation. Pre- 
viously, the two-loop renormalization of the pure Yang-Mills theory in the background 
field method was first considered in || using the MS scheme. Fermionic contributions 
were added later in [K|. Calculations in the background formalism for an arbitrary 



value of the gauge parameter were presented in [11|. In the standard approach the eval- 



uation of the mass dependent QCD /5-function at the two-loop level was performed 
in ||12||. In the latter publication only an approximate expression for the two-loop coef- 



ficient was given. Because of the complexity of such calculations the result of [12|] was 
not confirmed by any other group until now. The background field method provides 
the easiest way to calculate a mass dependent /3-function because here only propaga- 
tor diagrams need to be evaluated. A general method for the evaluation of two-loop 



propagator type diagrams with arbitrary masses was recently proposed in ]TB|,|T^]. 
Our paper is organized as follows: in Sec. 2 we describe the calculation of the back- 
ground field propagator, from which we obtain the RG and the effective running cou- 
pling in the BF — MOM scheme in Sec. 3. The relationship between the BF — MOM 
and the MS coupling is presented in Sec. 4 (analytical) and Sec. 5 (numerical). For 
completeness we include a formula for the bare BF propagator in Appendix A. The BF 
Feynman rules and the BF propagator diagrams are included in Appendices B and C, 
respectively. 



2 The background field propagator 

To regularize divergences we will use the dimensional regularization procedure in 
d = 4 — 2e dimensions. In the background field approach, we only need to calcu- 
late the background field renormalization constant Za in order to obtain the charge 
renormalization constant. The complete list of two- loop diagrams as well as the Feyn- 



man rules in the background field approach may be found in || - |TjJ . Za is determined 
by renormalizing the background field propagator according to 

l + n(QV,{m?}) l+n„(Q 2 ,^,{m?}) ' 1 ' 

where Q 2 = —q 2 and /i is the subtraction point. Bare quantities carry an subscript - 
In the MOM scheme the condition 

U(Q 2 , f , 2 ,{m 2 })\ Q2= ^ = (2) 

is imposed on the renormalized self-energy function. The renormalized mass rrii in 
our calculations is defined as a pole of the quark propagator. In the MOM scheme Za 
and therefore also the RG /3-function depend on the gauge parameter £. The gauge 
parameter is renormalized by 

6> = ^3, (3) 

where Z 3 is the renormalization constant of the quantum gluon field. To circumvent 
problems connected with the renormalization of the gauge parameter we have chosen 
the Landau gauge £ = 0. 

We repeated the calculation of all two-loop diagrams in the background formalism 
keeping non-vanishing quark masses. All calculations have been performed with the 



help of FORM [15j using the algorithm described in |L3| . Our results agree with those 
presented in [||,[10],[n|] for the limit of massless quarks. The sum of all unrenormalized 
diagrams for an arbitrary value of the gauge parameter is given in the Appendix. 
The renormalized self-energy amplitude Tl(Q 2 ) has the form: 

n(Q 2 ) = (^)u 1 + m 2 Ui + - (4) 



where 



An J \AttJ 



Ui(Q 2 /^ 2 , {m 2 /^ 2 }) = H Ca In ^ + T F £ ( IF ( ^ ) - IT ( 
U 2 (Q 2 /^ 2 , {m 2 /^ 2 }) = In ^ + T F £ ( U 2 ( ^ ] - U 2 ( ] ) . (5) 




As usual, Ca, Cp and T F are the group coefficients of the gauge group and n F is the 
number of flavors. 



The results of our calculations for 111 2 read 



n 2 (Qpj = ^±M.[{c A + AC F ) a{z) - {C A - 2C F )(1 - 2z) I(z)) 

+ — {39 + 31^ (z) - [Az 2 + 13Az + 57 - 12(2 - hz)zG(z)} (1 - z)G(z) 
Qz ^ 

+ 2[z 2 + I8z + 9 - 3(3 + 8z)(l - z)G(z)} ln(-4z)} C A 
+ |- {13 - [6(3 + 2z) + (7 + 8z- A8z 2 )G(z)](l - z)G(z)} C F , (7) 



where 



Q 1 -<, 



4m 2 J 1 -l/z+1 

denote the kinematic variables and 



y - 1 



I(z) = 6[C 3 + 4 Li 3 (-2/) + 2 Lia(j/)] - 8[2 U 2 (-y) + Li 2 (y)] my 

-2[21n(l + y)+ln(l-y)]ln 2 y, (10) 
= 6C3 - 6 Li 3 (y) + 6 hxy Li 2 (y) + 2 ln(l - y) In 2 y , (11) 

<r(z) = 3—^- {2 Li 2 (-y) + Li 2 (y) + [ln(l - y) + 2 ln(l + y) - ^ In y] lny}(12) 

are our basic integrals. The functions ^3 (^) are master integrals considered 

in 0,|l3. 

Setting Ca — 0,C F — T F — l,n F — 1 and taking the limit /1 2 — ► we reproduce the 
well known result for the photon propagator [P7|,|18"| in the on-shell scheme. 
At large Euclidean momentum Q 2 = —q 2 we find the asymptotic forms 

/Q 2 W-oo 4, Q 2 (m 2 \ (m 2 \\ Q 2 



n 2 ( ^) Q -|(5CU + SC F ) ^$ + l(C A + S6(C a - 2C F )( 3 ) 

6(3CU-8C7 F )(^)ln£L + .... (13) 



With these results at hand we are able now to obtain the mass-dependent two-loop 
/3-f unction. 



3 The RG equation and the effective coupling 

In the BFM the RG /^-function is given by 

d ( a s \ ,. d , „ „ /a s \ 2 „ ' 



]i m a s ^\nZ A = -P (^-) (14) 



2 

^ V V4tt; -"-dp—" " u V4vry " x V4vry 

and hence the coefficients of the /3-function may be simply obtained by differentiating 
(H) and ([?]). The results read 

A = ^- r F f>(^), (15) 



1=1 



where 



&d(4) =l + f(l-G(x)) 
m / 2a; 



6, (£) = [16(1 - + (1 + 8^) C . 4 ]^L - ±(C A - 2C,)/<*) 

2 4 
+ — /i 4) C A + [(1 + 3a; - 10s 2 + 12a; 3 )C A - 3(3 - 3s - 4a; 2 + 8x 3 )C F ]—G 2 (x) 
3x 6X 

- [(147 -Ax- 100a; 2 + 8x 3 )C A + 168(1 - x)C F + 6(9 + Ax) ln(-4x)C A ] — G(x) 

9x 

+ [(99 + 62x)C A + 12(11 + 3x)C F + 2(27 + 2Ax - 2x 2 ) h\(-Ax)C A \ — , (16) 

9x 

with a; = — /i 2 /(4m 2 ). This is our main result. 



In Ref. |L2| the /3-function for QCD ( C A = 3, Cp = 4/3, T F = 1/2) was evaluated in 
the standard approach with a renormalized coupling constant defined via the gluon- 
ghost-ghost vertex in the Landau gauge taken at the symmetric Euclidean point. The 
authors presented only an approximate result for the function 

MC\ - 3ft 

Bl(r) = + < 17 » 

which corresponds to our function b\(r) and which they parametrized as 

(-0.45577 + 0.26995r)r 



BAr) 



1 + 2.1742r + 0.26995r 2 



with r = /i 2 /m 2 . As asserted in |12 the parametrization flT8|) has the maximum 
deviation from the true value in the entire range < r < oo smaller than 0.005. We 
find that the difference between our expression and ( |I8"D in the same region is less than 
0.015, which is also very small. This is somewhat surprising, since we are comparing 
couplings in different schemes. 



For a mass-dependent renormalization schemes the RG equations 



in general can be solved only by numerical integration. However, an approximate solu- 
tion for the mass dependent effective QCD coupling was proposed in ||LS| ,|2~0[j. Indeed, 
at the two-loop level the expression 

Qs{Q2) = 1 + a s /(47r)C/ 1 + a fl /(47r)(C/ 2 /C/ 1 ) ln(l + a s /(4vr)f/ 1 ) ' (20) 

with U\ t 2 given in (|5|), correctly sums up all leading as well as "next-to-leading" terms 
a s ll2{<y s Ui) n though it is not an exact solution of the two-loop differential RG equation. 
We will compare (|20|) with the result of the numerical integration of the RG equation 
below. 



4 BF — MOM coupling in terms of the MS coupling 

Let us define the auxiliary functions 

„ , , 20 4, 

z X i = -Hi {n) - — - - k , 

z 2i = -n 2 (r,) - (f + jk) Ca - (f + 4/,) C F , (21) 

where 

Zj In /"j , 7"j „ . 

mf 

For later use we note that for light fermions, utilizing the expansion (|13|), we obtain 
z u = -y + <3(m 2 //i 2 ) , 

** = - + 8 Cs) ^ - (f - 16 Ca) C F + OK/// 2 ) . (22) 

The relationship between the renormalized coupling constants may then be written in 
the form 

h = ^^ = H(h,^)=h + k 1 {^)h 3 + (h 2 {^) + ^))h s + --' (23) 



where 



205 np 

k ^ 2 ) = ~na^" A + T f H z ^ ' 
36 i=i 



H^=( 2 ^-fCs)cl + T F Y i ^, (24) 



and 



h=h _ a s MS 



Differentiating the relation with respect to /i we obtain: 

2 dh - dH(h,fi 2 ) 2 dH(h,n 2 ) 

a* ^ = /wo = —Q^—m + „ 

= — Po mom h 2 — Pi mom h 3 — ■ ■ ■ , (25) 
where /3(h) is the (3 function in the MS scheme: 

(3(h) = -p h 2 - Pih 3 , (26) 

with 

(3q = ~2~Ca — -^T F n F , 
34 20 

A = y Ci - —CaTfUf - AC F T F n F . (27) 
From the above equation we obtain 



_9 

9// 

9// 

As it should be, in the massless limit the /3-functions agree. 



A) MOM — A) — /i 2 ^^-fci(/i 2 ), 

A mom = A - /J?-^k 2 (iJ?). (28) 



5 BF — MOM versus MS coupling: numerical aspects 



For numerical studies we use the following pole quark masses [23 



m u ~ m d ~ m s ~ ; m c = 1.55GeV ; = 4.70GeV ; m t = 173.80GeV . 

(5) 

For the strong interaction coupling we take cr^ = 0.12±0.003 at scale M^=91.19 



GeV p4 |. In Fig. 1 we show that Shirkov's formula (|20| ) provides an excellent ap- 
proximation to the exact solution of the two-loop RG equation. At sufficiently large 
scales the mass effects in the /3-function are small and we expect no large numerical 
differences between different schemes. This is illustrated in Fig. 2, where the evolution 
of the running couplings is shown for a common start value of a s = 0.12 at the scale 
Mz = 91.19 GeV. Only the space-like E = y/—q 2 is considered. We see that the mass 
effects are of comparable size as the 3-loop contribution pll, EI3] in the MS scheme 
(see Tab. 1 given below). The MS results were obtained by adopting the Bernreuther- 
Wetzel (BW) [E2J matching scheme between the effective theories with different flavors. 
We checked that utilizing Marciano (M) matching ||, instead, leads to answers some- 
what closer to our MOM results. Only the latter one exhibit the correct physical mass 
behavior. 
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Figure 1: Evolution of a s in the BF — MOM scheme normalized to a s = 0.12 at the scale 
Mz = 91.19 GeV. The dotted line represents the approximation by Shirkov's formula. 
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Figure 2: Comparison of the a s evolution in the space-like region normalized to a common 
value a s = 0.12 at scale Mz = 91.19 GeV. The dotted, dashed, dash-dot and the dash-dot- 
dot-dot curves show, respectively, the one-loop, two-loop, three-loop and the four-loop 
MS evolution for BW-matching. The full line represents the exact BF — MOM running 
coupling. 



Although BW-matching seems to be better justified from a field theoretical point of 
view, it leads to "threshold jumps" which of course are not physical in the space-like 
region. In contrast M-matching assumes continuity of a s across the matching scale 
("thresholds"). 

While there are no really large numerical differences in the /3-functions, i.e., the deriva- 
tives of a s with respect to fi, down to moderately low scales, there are large y,- and 



hence mass-independent terms in the relationship between the coupling constants (23Q, 
as follows from (123) and (E5J), and as we can see in Fig. 3. A large constant shift in 



55 



1.00 



2.72 



7.39 



20.09 



E (GeV; 



8.0 

7.0 

6.0 

5.0 

4.0 

3.0 

2.0 

1 . OH 

0.0 



2-loop MS M 

2-loop BF-MOM normalized to MS at M z 



2-loop BF-MOM 



2-loop MS M via BF-MOM ,\\{ 

..„„,iill#^ 




a s (M z )=0. 120 ± 0.003 (LEP 98) 
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Figure 3: Comparison of a s in the BF — MOM and the MS schemes with cr = 0.12 at 



scale Mz=91.19 GeV and a s bf-mom at this scale calculated using (23). The dotted line 
is the MS coupling calculated from a s bf-mom^) by inverting (|23"|). 



a s of about plus 14% at Mz is obtained when we go from the MS to the MOM scheme. 
In principle, this does not affect the prediction of physical observables. However, the 
scheme dependence which is due to truncation errors of the perturbation expansion is 
different for different renormalization schemes. The shaded area of Fig. 3 reflects the 
theoretical uncertainty at the two-loop level which shows up in the comparison of the 
two schemes. Below about 1.15(2.92) GeV the one-loop correction ki(y 2 ) h in (f2~3|) 
exceeds the leading trivial term by 100(50)% and the perturbation expansion cannot 
be applied any longer (see Fig. 3). 



The occurrence of the disturbing large numerical constants in the relationship between 
the renormalized couplings belonging to different renormalization schemes is not a 



peculiar feature in the relation between the MOM and MS schemes. Similar worrying 
lare terms, long time ago, were the reason for replacing the original MS by the MS 
scheme Q, which are related by a simple rescaling of the scale parameter fi. Other, 
more sophisticated, examples of eliminating leading terms by rescaling were proposed 
in Ref. [p5fl . Also, for the comparison of non-perturbative calculations of running 
couplings in lattice QCD with perturbative results, the adequate choice of a relative 
scale factor turns out to be crucial ||26|| . The rescaling usually leads to dramatically 
improved agreement. A condition for the rescaling to make sense is that the /3-functions 
of the two schemes under consideration do not differ too much numerically. For our 
two schemes this condition is fairly well satisfied (see Fig. 2). In fact at higher energies 
the /3-functions become identical. Such a rescaling procedure thus looks natural if we 
tune the running couplings to agree with good accuracy at high energies. This can be 



achieved as follows: While ( p3|) reads {k 2 = k 2 + k\) 

hi/j 2 ) = h(/2 2 ) + hh 2 ^ 2 ) + ~k 2 h 3 (i2 2 ) + 0(h 4 ) (29) 
we may absorb the disturbing large term ki into a rescaling of \i by a factor x such 



that 26 



h{{x^f) = h(fi 2 ) + + 0(h 3 ) . (30) 

Expanding the RG solution fl2"U| ) we have (U 2 = U 2 — Uf) 

h((x fi) 2 ) = h(fi 2 ) - EM4 K/VI^V) - U 2 (x 2 , {m 2 J^})h 3 (ji 2 ) + 0(/i 4 ) 

= h(f, 2 ) + (h - Ut(4, {m 2 /^ 2 })) h 2 (f, 2 ) (31) 

+ (k 2 - 2U t (x 2 , {my^h - U 2 (x 2 , {m 2 /^ 2 })) h\^ 2 ) + O(h') 

and the rescaling factor xq is determined by the equation 

k^U^xlim 2 /^ 2 }). (32) 

In our mass dependent scheme we require this to be true only at very large scales 
f-i 2 ^> m 2 for all flavors / including the top quark. This convention is simple and most 
importantly, it does not conflict with the manifest decoupling property of the MOM 
scheme. As a consequence we obtain a running coupling which depends very little on 
the scheme at large energies, a property which looks most natural in an asymptotically 
free theory like QCD. For the BF — MOM scheme the rescaling factor xq is determined 
by 

HA) = (i^Ca - y TfUf ) I {y Ca ~ t TFnp ) = 125/84 (33) 

for QCD with rip = 6 flavors. Numerically we find xq ~ 2.0144. 



In order to check whether the above rescaling makes sense, we must inspect the change 
of the 2-loop coefficient in the rescaled relationship ( |3"2"D between MOM and MS . 
Indeed, the rescaling changes the coefficients from k\ ~ 10.42, k 2 ~ 126.35 to 
k\cE = 0, k 2e ff ~ —32.46 and thus we get a substantial improvement for the next 
to leading coefficient too, as it should be. We note that the rescaling improved MOM 



perturbation expansion at low energies does not any longer deviate substantially from 
the MS results. Of course, only the appropriate higher order calculations of observ- 
ables in the BF — MOM scheme could reveal the true convergence properties of the 
perturbation series in this scheme. 

In the MOM scheme the energy scale comes in by a momentum subtraction and the 
location of the thresholds of course cannot depend on the rescaling "reparametriza- 
tion". This means that actually the scale must be changed in the MS scheme, where 
the scale parameter \x enters in a purely formal way and "thresholds" are put in by 
hand for switching between the effective theories of different numbers of flavors. Since, 
conventionally, /x in the MS scheme has already been identified with the cm. energy, 
for example, in the LEP determination of a s (Mz) which we use as an input, we have 
to apply the rescaling to the MOM calculation. As the thresholds must stay at their 
"physical" location, i.e., Am 2 /q 2 must remain invariant, we have to perform the scaling 
simultaneously to the energy and the masses. 



The result from utilizing this rescaling procedure is displayed in Fig. 4. The large 
deviations seen in Fig. 3 have disappeared now. The sizes of effects are still illustrated 
by what we observe in Fig. 2 except that the initial values at Mz differ. In Fig. 4 we 
have recalculated the input values of af\Mz) as a function of the perturbative order, 
assuming the observable R(s) to have a given experimental value. R(s) is the ratio of 
hadronic to leptonic e + e~-annihilation cross sections at sufficiently large s, from which 
a precise determinations of a s (s) is possible. At our reference scale M z we may use 
perturbative QCD in the massless approximation P7[ 



R{s) = (l + a + Cl a 2 + c 2 a 3 + ■■■) (34) 

/ 

where Qf denotes the charge of the quark, a = Ah = a s (s)/ir, and 
ci = 1.9857- 0.1153 n F 

c 2 = -6.6368 - 1.2002 n F - 0.0052 n% - 1.2395 (^ Q/) 2 /(3 Q)) 



in the MS scheme, with n F = 5 active flavors. 



Some concluding remarks: We have investigated a MOM renormalization scheme in 
the background field gauge at the two-loop order in QCD and shown that a substantial 
scheme dependence is observed relative to the MS scheme, unless we apply a suitable 
rescaling. These findings are in accord with earlier investigations at the one-loop ||^8 
and two-loop 



12 level. 



Mass effects in any case are non-negligible at a level of pre- 
cision where also higher order corrections are relevant. The calculation in full QCD 
includes the exact mass effects and is smooth and analytic at all scales and in particular 
across thresholds. It thus avoids problems with the MS scheme addressed in a recent 
article by Brodsky et al. [2{J which were cured by an analytic extension of the MS 
renormalization scheme. 



We note that the use of the BF — MOM scheme, particularly when using the compact 
form obtained for Shirkov's approximation, is much easier in practice because decou- 



pling is manifest at any threshold and there are no matching conditions to be imposed. 
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Figure 4: Comparison of a s ((xo-E') ) in the BF — MOM and a s (E ) in the MS scheme with 

(5) 

input values = 0.120 at scale M^=91.19 GeV and a s bf-mom = 0.1189 obtained for 

the rescaled energy x M z ~190.90 GeV (x ~ 2.1044). 



We emphasize that the scheme and scale dependence of perturbativ QCD predictions 
is not a matter of the order of perturbation theory alone but may depend substantially 
on other details like the kind of matching condition applied in the mass independent 
MS schemes or the threshold and mass effects in MOM schemes. The following table 
(Tab. 1) may illustrate the kind of uncertainties we expect to encounter. We find that 

Table 1: Comparison of predicted a s values at the masses of the T, J/ip and r. Here 
we adopt a common input value a s (Mz) = 0.12 for the MS scheme independent of the 
perturbative order. 



scheme 


a s (M z ) (input) 


a s (M T ) 


a s (M J/v ,) 


ot s {M T ) 


MS 2-loop (BW) 


0.120 


0.179 


0.260 


0.354 


MS 3-loop (BW) 


0.120 


0.179 


0.262 


0.364 


MS 4-loop (BW) 


0.120 


0.179 


0.263 


0.368 


MS 2-loop (M) 


0.120 


0.179 


0.258 


0.348 


MS 2-loop via BF - MOM 


0.120 


0.168 


0.211 


0.254 


(BF - MOM 2-loop 


0.120 


0.180 


0.265 


0.358) 


BF - MOM 2-loop 


0.137 


0.222 


0.372 


0.605 


BF - MOM 2-loop rescaled 


0.121 


0.181 


0.260 


0.345 



the 2-loop MS value at the r-mass M T is a s = 0.254 when we switch from MS to 
BF - MOM at the Z-mass M z use BF - MOM evolution down to M T and switch back 
from BF — MOM to MS . Standard (direct) MS evolution depends on the matching 
scheme utilized (BW or M) and for BW(M)-matching yields a s = 0.354(0.348), such 
that via BF — MOM we get a value which is lower by 0.100(0.094). However, the 
BF — MOM value obtained with the rescaling is 0.345, not very different from its MS 
value. The Particle Data Group p3| quotes a s (M T ) = 0.35 ±0.03 for the experimental 



value obtained from r-decays (see also |j30[ ) 



Experience with many physical applications of the MS scheme somehow established 
this scheme as a preferred one, in the spirit that this prescription is better than others 
in the sense that it leads to reliable perturbative predictions for many physical observ- 
ables. In our opinion it remains unclear whether a prefered scheme exist. The problem 
is the appropriate choice of scale. We advocate here to take more serious the physical 
mass dependence. In order to get a better understanding of the scheme dependences 
we need more calculations in different schemes. 
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Appendix A: Bare BF propagator 
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d = (3d - 8)£ [-(d - l)(d 2 - 9d + 22)(d - 4) 2 £ 2 - (15d 5 - 256d 4 + 1685d 3 
-5292d 2 + 7744d - 3960)£ + 33d 5 - 542d 4 + 331 Id 3 - 9378d 2 + 12448d - 6608 
+ 27112d - 8128 - 33312d 2 - 51d 6 + 916d 5 + 20016d 3 - 6169d 4 , 



c 2 = (d- l)(d - 4) 2 [(d - 4)£ + 6(2d - 7)]£ 3 + 2(lld 4 - 144d 3 + 677d 2 - 1296d + 743)£ 2 
- (84d 4 - 854d 3 + 2994d 2 - 4304d + 2384)£ + 49d 4 - 403d 3 + 1106d 2 - 1392d + 664, 

/i = (d - 2) [(d 2 - 7d + 16) z - (d - 5)(d - 4)]z - 2, 



f 2 = (d- 2)[(d - 4)f - 7d + 12]z + (d - 4)£ 2 + 2((d - 2)z + 4)(3d - 10)£ - 7d + 16, 
h = 2(d - 2)(d 2 - 5d + 8)z 2 + (d - l)(d - 3)(d - 4) 2 z + d 3 - 6d 2 + 5d + 8, 
U = (d-2)(l-2z)[(d-2)z + l], 

f 5 = 4[(d - 1)£ + 3d - 7](d - 4)z 2 - 10[(d - l)(3d - 8)£ - 7d 2 + 23d - 28]z 
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- (d - 7)(9d - 41) (2d - 7) (3d - 8)£ + (3d - 8)(26d 3 - 275d 2 + 892d - 787), 

f 9 = 2(d - 5)(d - 3)(d - 4) 2 [(d - l)f + 3d - 7}z 3 

- (d - l)(d - 4)(19d 3 - 208d 2 + 711d - 754)£ - 6266d + 345d 4 + 2888 - 23d 5 
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/ 10 = (d - 2) [(d 2 - 5d + 8)z + d 3 - 7d 2 + 16d - 14]. 
In the above formulae we have used the following notation: 
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All parameters are the bare one's, Zi = q 2 /(Am 2 ) and the coefficient functions f n = 
fn(zi) are functions of Z; L . 



Appendix B: The BF Feynman rules. 

In addition to the conventional QCD Feynman rules we have: 
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All momenta are taken to be outgoing. 
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